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Abstract 

In the preceding paper, the structure and thermodynamics of a given quantum system was 
represented by a corresponding classical system having an effective temperature, local chemical 
potential, and pair potential. Here, that formal correspondence is implemented approximately 
for applications to two quantum systems. The first is the electron gas (jellium) over a range of 
temperatures and densities. The second is an investigation of quantum effects on shell structure 
for charges confined by a harmonic potential. 



I. INTRODUCTION 



In a companion paper pQ a method was described that would allow application of strong 
coupling classical many body methods to calculate properties of equilibrium quantum sys- 
tems. Within the grand ensemble for equilibrium statistical mechanics, the thermodynamics 
and structure are obtained as functions of the temperature T, the local chemical potential 
\i (r) =ii — 4>ext (r), and a pair potential 0(r,r'), (where (p ext (r) is an external single par- 
ticle potential). A grand ensemble for a corresponding classical system is characterized by 
an effective temperature T c , an effective local chemical potential \i c (r) = fi c — <fi c ,ext( r ), 
and an effective pair potential c (r,r'). These three classical parameters are fixed by three 
correspondence conditions: equivalence of classical and quantum pressures, densities, and 
pair correlation functions. An approximate inversion of these formal definitions to obtain 
T c , /U c (r), and c (r, r') was described within classical liquid state theory. The objective 
here is to illustrate this approach for two applications. The first is to calculate the pair 
correlation functions for the electron gas (jellium), the prototypical test bed for quantum 
correlations [21 E] , over a wide range of temperatures and densities. Corresponding ther- 
modynamic properties can then be calculated in terms of these correlation functions. The 
second application is to harmonically bound charges in a trap [I]. Specifically, the role of 
quantum diffraction and exchange as a mechanism for shell formation is investigated. While 
both systems have been studied extensively at both very low (ground state) and very high 
(plasma) temperatures, the relevance here is a method that applies across the intermediate 
domain. 

The use of effective pair potentials to include some quantum effects in classical methods 
like molecular dynamics simulation has a long history [5 J . A new phenomenological approach 
proposed by Perrot and Dharma-wardana [6] more recently goes a step farther to introduce 
an effective classical temperature as well. Applications of this extended approach to a variety 
of systems and properties over the past decade have met with remarkable success [7]. The 
present work can be considered as a parameter free formalization of this earlier work, and 
comparisons are discussed critically here as well. 

In the next section, the approach of reference [1] is applied to jellium. The effective 
classical pair potential is discussed and illustrated, and then the pair correlation function is 
calculated using the classical strong coupling hypernetted chain integral equations (HNC) [H] . 
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The dimensionless temperature relative to the Fermi temperature, t = T/Tp, is considered 
in the range < t < 10. The dimensionless measure of the density is r s = r /ciB, the mean 
distance between particles defined by Airr^/S = 1/n relative to the Bohr radius ag. It is 
considered in the range 0.5 < r s < 10. The results are compared to the above mentioned 
model of Perrot and Dhama-wardana (PDW), and to the extension of the self-consistent 
generalization of the random phase approximation, STLS [9], to finite temperatures by 
Tanaka and Ichimaru [10]. Comparison with limited diffusion Monte Carlo simulation data 
for t = [TT] shows good agreement. No benchmark simulation data at finite temperatures 
is available as yet. 



A second application is described briefly in Section |III[ to charges in a harmonic trap 
where classical strong Coulomb correlations produce shell structure pi] . These results are of 
interest for laser cooled ions in traps, where the quantum effects are expected to be weak, and 
for electrons in quantum dots at low temperatures, where the quantum effects are strong. 
In the classical case, strong Coulomb correlations are required for shell structure - they are 
absent in a mean field theory at any value of the coupling constant. The objective here is to 
explore possible new origins for shell structure due to quantum effects. It is shown that the 
quantum mean field theory, without any Coulomb correlations, leads to shell structure due 
to diffraction effects modifying the Coulomb interactions and/or exchange effects modifying 
the effective trap potential. 

The results are summarized and discussed in the last section, where a modified version 
of the simple PDW model is proposed. 



II. APPLICATION TO UNIFORM ELECTRON GAS 



The interacting electron gas is a one component system of charges with Coulomb in- 
teractions embedded in a uniform neutralizing back ground. The uniform electron gas, or 
"jellium", provides an important model system to discuss correlations and quantum effects 
in real physical metals, solids, and plasmas |2J. The classical limit is known as the one 
component plasma. There are two parts to this section. First, the parameters /3 C = 1//cbT c , 
/x c (r), and c (r, r') for the effective classical system are determined approximately from their 
definitions in reference [1] (for the uniform jellium /i c (r) — > fi c and c (r, r') — > c (|r — r'|)). 
Second the resulting classical system is applied to calculate the structure and thermody- 
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namics of jellium from the classical HNC liquid state theory. 

To summarize the results of pQ, the classical parameters are defined as following. The 
classical temperature is obtained from the correspondence condition of equal pressures and 
the classical virial equation 



Pc _ PcPc _ n_ 
13 ~ p P ~ p P 



ft 

1 - - / drh(r)r ■ V/3 c c (r) 



( . , - . - • , . (1) 

The replacement of the pair correlation function g(r) by the hole function h(r) = g(r) — 1 
occurs because of the uniform neutralizing background. The classical activity /3 c /x c is given 

by nsi 

(3 cf i c = In (n c Xl) -n J dv (c{r) + /3 c c (r) - l -h{r) (h(r) - c(r))J , (2) 

where c(r) is the direct correlation function defined in terms of h(r) by the Ornstein-Zernicke 
equation [S] 

c(r) = h(r)—n J dr'c(\r — r'\)h(r') . (3) 
Finally, the pair potential is obtained from the inversion of the HNC equation 

/3 c c (r) = - In (1 + h(r)) +h{r)-c (r) . (4) 

The classical pair correlation functions on the right sides of these expressions have been 
replaced by the quantum functions, according to the third correspondence condition. Hence 
these classical parameters are determined by quantum input. 

The practical approach is to provide the essential quantum input by specifying h(r) in 
some approximation. Equations Q and ^ then determine (3 c (f) c {r) and c(r), and with 
these known (3 C / (3 and /3 c /i c can be calculated. The objective here is to propose a simple 
approximation for practical application. 



1. Classical potential f3 c (j) c {r) 

The dominant exchange effects are already present in the ideal gas calculation described 
in reference [T|. Therefore it is convenient to write /3 c c (r) in the form 

/3 c c (r) = (0 c (r)) (o) + A(r), (5) 

where (0 c (r)) < -°' ) is the ideal gas Pauli potential and A(r) denotes the contribution to the 
effective potential from the Coulomb interactions. In the classical limit A(r) — > (3q 2 jr. 
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Another exact limit is the weak coupling limit for which the direct correlation function 
becomes proportional to the potential, or stated inversely, 

PMr) -> -c(r), (/3 c c (r)) (o) -> -c^\r). (6) 

Thus a possible approximation incorporating this limit is 

/3 c c (r) (/3 c c (r)) (O) - (c(r) - c( )(r)) W , (7) 

where (c(r) — c^(r))^ w denotes a weak coupling calculation of the direct correlation func- 
tions from the Ornstein - Zernicke equation (|3]). For the classical OCP (Coulomb potential) 
this yields the Debye - Huckel approximation to h (r). Here it is required that this should 
yield its quantum counter part, the random phase approximation (RPA) [3]. The weak 
coupling calculation from the Ornstein-Zernicke equation is then 

c ( r )H = h RPA ^_ n J dr > (c( | r _ r /|))W h RPA^ (g) 
This has the solution 

c(r )H = I [ * e -*-rQ|-i. (9) 



n 



(2tt) 3 S RPA {k) ■ 



Here S RP (k) is the RPA static structure factor 

S RPA (ty = i + n f dre ik-r h RPA( r y ( 1Q ) 



Finally the modified Coulomb potential A (r) in ^ becomes 

AW ^/^' kr (^-s4))- (u) 

The definition of S RPA (k) in terms of the RPA dielectric is given in Appendix [Aj 

Several limits of A (r) = A (t, r s , r*) are established in Appendix |A} For large r* = r/r 
it behaves as r* _1 

UmA(t,r s ,r*) ^T e (t,r s )r*-\ (12) 
where T e (t, r s ) is an effective Coulomb coupling constant 

Te (t,rs) = /3^coth(/3na; p /2) r ' T = V (13) 
Here a; p = \J Airnq 2 /m is the plasma frequency. The dimensionless parameter is /3hoo p = 
(4/3) (2-\/3/vr 2 ) 1 ^ 3 ^Jr~s/t so for fixed r s the high and low temperature limits are 

r -"V VT 1 , (14) 

(§r a ) , /3^ p » 1 



This asymptotic Coulomb form is exact and its coefficient follows from the fact that the RPA 
incorporates the exact perfect screening sum rule [12]. It is illustrated for r*A (t,r s ,r*) in 
Figure [T] at r s = 5 for several values of t. 

Also shown in this figure are the results from the PDW classical potential 

WMr)f DW = WMr)) {0) + A PDW (r) . (15) 

The Pauli potential (fic4>c(y))^ as in is the same as in (JsJ) but its correction A PDW (r) is 
given by the Deutsch regularized Coulomb potential [13J 

„2 , . / oPDW fc2\ V 2 x 
A PDW (r) = pPDW<L f 1 _ e -r/A™\ X PDW = I P ^ \ = A 

Here /3 PDW = l/k B T PDW , where T PDW = (T 2 + T 2 ) 1/2 . The only free parameter, 
To = To (r s ), is fit by requiring that the classical correlation energy matches the quantum 
exchange/correlation energy obtained from quantum simulation at T = 0. The fit given in 
reference [H] is 

a + x /r s + cr s 



with a = 1.594, b = —0.3160, and c = 0.0240. It is seen that the PDW model is quite similar 
to the approximation defined here at r s = 5. Greater discrepancies occur for both larger 
and smaller r s except at higher temperatures. Further comments on this comparison are 
given below. 

For r* << 1, A (t, r s ,r*) approaches a finite value 

1 f f 1 1 \ 

A(Vr - 0) = ^(2^(s^-^w)' (18) 

The integral converges because the static structure factors for large k approach 1 as k~ 4 

due to quantum effects (cusp condition [15]). The Coulomb singularity is therefore removed 

in the effective classical pair potential. Finally, another limit obtained in the Appendix [A] 

is that for large r s and large t (low density, high temperature) in which case the Kelbg 

potential [2] is recovered 

lim A(t,r s ,r*) ->A K (t,r*) = — (l - exp(27r ^r* 2 ) + V2n^r*(l - erf (V^^r*))] , 
t,r s »l r* \ X K Xk Xk J 

(19) 

with Afc = A/ y/2n. The Kelbg potential is the exact weak coupling effective classical potential 
determined from the two particle electron - electron density matrix [14J. This limit is 
approached to within 10 percent at t — 10 and 1 < r s < 10. 




FIG. 1: Demonstration of crossover for r* A(t,r s ,r*) to Coulomb with effective coupling constant 



T e (t,r s ) given by eq.(13), for r s = 5 and t = 0.5, 1, 10. Also shown are the corresponding results 
for r*A PDW (t,r Sl r*). 

2. Classical effective temperature, chemical potential 

The approximate temperature and chemical potential equations are obtained in a similar 

way 

(3 C = $0) + (pRPA _ flf*.<») . (20) 

Pdic = (/W 0) + (M RPA - (M RPAM , (21) 

where (3^ and (f3 c /j, c )^ denote the ideal gas results of p], and from |IJ and (2) 

"l - f / dxh RPA {r)v ■ V (PMr)) RPA] 



n 

RPA 



pRPA 



(22) 



(M RPA = I In (^y) + ln ( nXi ) RPA + l n J drh RPA (r) (h RPA (r) + (/3 c c (r))^ A ) . 

(23) 

The RPA results for p RPA and (n\ 3 ) RPA are computed from the Pade fits of reference [2]. 

A peculiarity of jellium is the possibility for the pressure to become negative at large r s 
and small t, conditions for which the equal pressures correspondence condition cannot be 
imposed. For real systems, the pressure is positive as follows from the convexity of the free 
energy as a function of the volume. This convexity does not hold for jellium [16J. To be 
more explicit, it is first noted that the pressure is an increasing function of t so its minimum 
value occurs at t — 0. Figure M shows the quantum p RPA at t = as a function on r s . It 




-0.03 J 



FIG. 2: Quantum RPA pressure prpa at i = as a function of r s . 



is seen that p RPA (t = 0) vanishes for r s ~ 4, and becomes negative for larger r s . Thus for 
r s > 4 the pressure p RPA (t) vanishes at some temperature t (r s ). Then from (20) and (22), 
the effective classical temperature vanishes at t Q 

" RPA (to) 



T c (t (r s )) 



n 



1 - | / rfr/i RPA (r)r ■ V (f3 c <p c {r)) 



RPA 



0. 



(24) 



For t < to(r s ), p (t) < 0. However, it is found that the denominator of (24) remains 



positive. Since the classical temperature must be positive this indicates that the equivalence 
condition, p c \ = p, can no longer be realized. Therefore for jellium, this one of the equivalence 
conditions should be replaced by a different condition (e.g., equivalence of internal energies). 
Instead, the analysis here is restricted to t > t (r s ) to assure positive pressure. Figure [3 



shows t c = T C /T F as a function of t calculated from (20) for r s = 0,1,3,4, and 5. Figure [4 



shows the corresponding results for [i c /E F calculated from (21) 



A. Radial distribution function and thermodynamics 

With the parameters of the effective classical system determined approximately above 
they can be used in an accurate classical many-body method to combine the quantum 
properties of these parameters with classical strong coupling effects (e.g., molecular dynamics 
simulation). This is illustrated here by using the full HNC integral equations Q and (j3]) 
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FIG. 3: Classical reduced temperature T c /Tp as a function of t for r s = 0, 1,3,4, and 5. 




FIG. 4: Dimensionless classical chemical potential fj, c /Ep as a function of t for r s = 1, 3, 4 and 5. 

specialized to calculate g{r) for jellium 

In g(r) = — /3 c c (r) + h (r) — c (r) , c(r)=/i(r)— n J dr'c(\r — r'\)h(r') . (25) 

Of course, these equations are also those used to define (3 c (j) c (r) so the analysis would seem 
to be circular. However, the approach has been to use an approximation to the HNC 
equations to determine /3 C C (here the weak coupling RPA limit) and then to "bootstrap" 
this information to solve the full HNC equations for a g(r) that goes beyond the input 
g RPA {r) to include classical strong coupling. One manifest improvement obtained in this 
way is positivity of g(r), already noted in |J|. In contrast, g RPA {r) becomes negative for 
small r at sufficiently large r s . 
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The determination of g(r) from (25) is straightforward using the method described in 
reference [T7j. Note that these equations do not use the equal pressure condition nor the 
value of (3 C . Hence they do not have the restriction to positive pressures and the associated 
restriction on r s . The results are shown in Figure [5] for the case of r s = 5 at t = 0.5, 1, 



and 10. Also shown are the results from the PDW map defined by (15) and (16). The two 
effective classical system predictions are quite similar (differences of the order of 10% occur 
for r* < 1). This is somewhat accidental to the value r s = 5, and larger differences are 
observed at both smaller and larger r s except at high temperatures. 




FIG. 5: Radial distribution function g(r*) for r s = 5 at t = 0.5, 1, 10. Also shown are the results 
of PDW. 




FIG. 6: Radial distribution function g(r*) for r s = 5 at t = 0.5, 1, 10. Also shown are the results 
of TI. 
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Other theoretical models for g(r) are based on the same dielectric formalism of the RPA, 
but including "local field" corrections. One of the earliest was the self-consistent STLS model 
[5] at T = 0, later generalized to finite temperature T by Tanaka and Ichimaru (TI) [TUj . 
Figure [6] shows the same conditions as Figure [5] for comparison with TI. The discrepancies 
are largest at lower t and most noticeably at small distances where TI becomes negative. 
Not shown is the RPA which is even significantly more negative in this range. Both RPA 
and its improved TI overestimate the size of the electron correlation hole [3] at larger values 
for 7%. 




r s =1 
r B =5 
r B =10 

PDW [r 5 =1] 
PDW [r s =5] 
PDW [r s =10] 
DMC [r s =1] 
DMC [r s =5] 
DMC [r s =10] 



FIG. 7: Radial distribution function g(r*) for t = at r s = 1, 5, 10. Also shown are results from 
PDW and diffusion Monte Carlo. 

The PDW g(r) is in excellent agreement with diffusion MC data at t — [6] even at 
large r s . Although the single parameter T is fixed by fitting the t = exchange/correlation 
energy from MC data it is nevertheless impressive that this provides good results for g{r) 
across a range both r s and r. Figure [7] shows a comparison of the results of the present 
analysis with those PDW and MC predictions at t — for r s = 1,5, and 10 [11]. The good 
agreement is quite surprising since there is no MC parametrization in the present analysis 
and all quantum input is via the RPA and ideal gas exchange. However, it is recalled that 
the RPA preserves the exact quantum mechanics of the perfect screening sum rule that 
governs the cross over to exact large r* Coulomb limit. This is discussed further in the last 
section. 
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1. Thermodynamics 



The predicted pressure, p c in atomic units, for the effective classical system is obtained 
from 



PcPc 



n. 



1 — —n I dxh(r)v ■ V/3 c c (r). 



(26) 



and the effective temperature (22). Figure [8] shows this as a function of t for r s = 1, 3, 4 and 
5. Also shown are the corresponding results for modified RPA (using the fits from reference 

TO- 




FIG. 8: Dimensionless classical pressure p c /(nEp) as a function of t for r s = 1, 3, 4, 5. Also shown 
are the corresponding modified RPA results. 



III. APPLICATION TO CHARGES IN A HARMONIC TRAP 

As a final application here consider N charges localized within a harmonic trap. The 
Hamiltonian is 

*-^ = £ir+2£ k^-rl - / ^(r)n(r), (27) 

with the local chemical potential given explicitly as 

yu(r) = n — -muj 2 r 2 . (28) 

The constant \x determines the average number of charges N. As a consequence of the 
harmonic potential the equilibrium average density profile for the charges is non-uniform 
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and depends only on the radial coordinate 

oo „ 

n ^J2 N / dr 2-dr N (n-rj^l e -^ H -^ N) \r x ..r N ) , (29) 



n(r) 

N=0 



where (rj..rjv| X |ri..rjy) is the N particle diagonal anti-symmetric matrix element in co- 
ordinate representation, and Q is the grand potential. The density profile in the clas- 
sical limit has been studied in detail, via simulation and theory In that case the 
dimensionless form depends on r s and t only through the Coulomb coupling constant 

1 /3 

T = f3q 2 /tq = (4/3) (2/37T 2 ) r s /t. For sufficiently large F the formation of shell struc- 
ture is observed in n(r). The objective now is to exploit this effective classical description to 
explore the effects of quantum diffraction and exchange via the proposed effective classical 
system. Only a preliminary investigation of new mechanisms for shell structure is described 
here, with a more complete discussion to be given elsewhere. 

The basis for the study is the HNC description for the inhomogeneous case, Eq. (37) of 
reference [1] 

In (n (r) A*) = c fJ, c {r) + J dr'c (2) (r, r | n)n (r) . (30) 

The classical studies of reference jl] made a further approximation to this expression, replac- 
ing the correlations for the inhomogeneous system c^(r, r" | n) by those for a corresponding 
uniform one component plasmas (OCP or classical jellium), c^(r, r" | n) — > c(|r — r"| , n). 
The results based on this approximation are found to be quite accurate except at very strong 
coupling. A partial theoretical basis has been given |19j . This approximation will be made 
here as well. 



It is convenient to rewrite (30) in a Boltzmann form with an effective potential U(r) 
defined by 



e -PcU(r) 



"< r > = J V«^T (31) 



so that (30) becomes 



N f 

0JJ(*) = ~Pc Mr) - ^) - dr , e _, cU{rl) J dv'e-^c{\v - r'| ,n). (32) 



Practical application of this result requires specification of the direct correlation function 
c(r,n) for jellium and the classical local chemical potential /z c (r) — \i c . The former is deter- 
mined from the equivalent classical calculation described in the previous section. The latter 
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is the effective classical trap potential corresponding to the actual quantum harmonic trap. 
It's approximate determination is described in the next subsection. 

The total number of particles appears explicitly. To introduce the density, it is necessary 
to assign a volume for the system. This can be defined as the volume of a sphere with radius 
Ro corresponding to a particle at the greatest distance from the center. At equilibrium the 
average density can be taken to be spherically symmetric so that the total average force on 
that particle is 

^ - mu 2 R = 0, R 3 = N-^—=. (33) 



This gives the average density to be 

3A 3mu 2 . A ^ 

^4^3 =4^T ( 34 ) 

In this way, the trap parameter mu 2 /q 2 is specified in terms of the density. 



A. Approximate form for fi c (r) 



Without quantum effects fi c — /x c (r) is just the harmonic potential of (28). Modifications 
for the effective classical form occur due to both diffraction and exchange effects. The 
exchange effects are dominated by those for the ideal Fermi gas in a harmonic trap. As 
a first approximation here, /i c (r) is replaced by that for an ideal gas of N Fermions in a 
harmonic trap. This is the inhomogeneous ideal gas considered in Section IV in pp. The 



local chemical potential is therefore obtained from (30) specialized to an ideal gas 

(/3 c /i c ) (0) (r) = ln(n (0) (r)A^) - J dr'c (0) (|r - r'|)n (0) (r')- (35) 
Then the effective potential U(r) becomes 

/W(r) = -ln(n<°>(r)A2)- J dr'c^(\r - r'\)n^(v')+^ c -N J dr'j£^^c(\r - r'\,n). 

(36) 

The direct correlation function for the uniform ideal Fermi gas, c^°\r, n), is known from the 
results of pQ. Furthermore, n™(r) has an explicit form in the local density approximation 
of Appendix B, reference [1] 

2 



»<"<r) = . (37) 
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where T is the same coupling constant as in eqn (13). The Fermi function / 3 / 2 (x) is given 
by 



dxx 2 z~ l e x + 1 



(38) 



The constant chemical potential fi^ is determined in terms of t, r s by the condition that 
the average number of particles is N 



N = 2 



roN3 
A 



dr*f 3/2 [e 



^( M (o)_i rr *a) 



(39) 



Hence n^(r)rg is also given by (37) in terms of t,r s . Figure |9j shows (3 c (J) cex t(r) = 
(/3 C // C )^(r) — (P c fj, c p°\0) obtained in this way for r s = 5 and t = 0.5,1,5, and 10. The 
classical trap potential is harmonic at large r, but there are significant deviations at the 
lower temperatures for 2 < r < 3. 




FIG. 9: Effective classical trap potential Af3 c /i c (r*)/T = (/3 C // C (r*) — f3 c fj, c (0)) /T as a function of 
r* for t = 0.5,1,5,10. 



B. Quantum effects on the mean field density profile 



The density profile can now be determined from (32) where the effective potential of (35) 
becomes 

P c U{r) = M c , ext {v) -N J dv' ^_* cU{rl/) c(\r - r'| ,n). (40) 
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(the potential c U(v) has been shifted by a constant to simplify the result). The effective 
classical external trap potential is 



/3 C 0c,ext(r) 



ln(n (0) (r)A 3 ) + J dv'c^(\v - v'\)n^(v') 



(41) 



Quantum effects result from the deviation of (3 c 4>c,ext{ r ) from the given harmonic potential, 
and the deviation of c(r, n) from its classical OCP form. In this section these two sources 
are isolated to explore the possibility of new origins of shell structure. To do so only the 



mean field limits of (40) and (41) are explored here. The mean field limit is defined by 
c(r, n) — > —f3 c (j) c (r) and S-°\r } n) — > — (/3 c c (r)) < -°' ) so (40) and (41) become 

e -PcU(r') 



P c U(v)^f3 c <p c , ext (v) + N / dv' 



(/3 c c (|r-r'|)) (O) + A(|r-r'|) 
dv' (/3 c c (|r-r1)) (O V°)(r')- 



(42) 



/3 c c ,e**(r) -ln(n(°)(r)A 3 ) -J dv' (/3 c c (|r - r'|)) w n^(v'). (43) 

The classical mean free limit corresponds to (/3 c (j)c( r ))^ — and A (r) = /3q 2 /r. There is 
no shell structure in this limit, even at very strong coupling. Instead, shell structure arises 
due to sufficiently large Coulomb coupling, such that c(r, n) differs from —(3q 2 /r inside the 
correlation length ro and the Coulomb singularity at r = is removed. This is shown in 



Figure 10 for T = 3. 




FIG. 10: Comparison of — c(r*) and Vk{t*) at t = 0.1,0.27 both corresponding to T = 3. Also 
shown is the Coulomb limit j3q 2 jr. 
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1. Diffraction effects 



To explore the effects of diffraction only in (|42j) and (43), the contributions from exchange 
are set to zero, i.e. {P c (p c {r))^ — > 0, /3 c C)ex . 4 (r) — >■ mu 2 r 2 /2, and A (r) — > A K (r) = Kelbg, 



eq. (19). Then (42) becomes 



WW -+ Pc\m^ + N J dv' j^^ pA K (|r - r'|) 



(44) 



This has the same form as the classical limit, except with the Coulomb potential replaced 
by the Kelbg form. The latter differs from Coulomb at short distances, for which it is finite 
at r = 0. Thus diffraction (without Coulomb correlations) leads to the same qualitative 
physical effects as classical Coulomb correlations. This is illustrated in Figure [10] where 
the Kelbg potential is evaluated at r s = 0.042, t = 0.1 and at r s = 0.11, t = 0.27 both 
corresponding to T — 3. For this reason it can be expected that the quantum diffraction 
mean field approximation can give rise to shell structure not present in the corresponding 



classical case. This is shown in Figure 11 for T = 3 and t = 0.1,0.5,1,2. A clear shell 
formation occurs at the two lowest temperatures, for which the diffraction regularization of 
Ak (0) is greatest. 



r=3 






— -Coulomb 




1=0.1 




1=0.5 




1=1 

— 1=2 











FIG. 11: Diffraction mean field approximate density profile for T = 3 and t = 0.1, 0.5, 1, 2. 
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2. Exchange effects 



Now return to the mean field form (42) and neglect all diffraction effects to study the 
effects of exchange only 



c U{r) -+ 0, 



c'rc,ext 



r) + N / dr' 



-PcU(T>) 



(^(Ir-r'^^ + ^lr-rr 1 ! , (45) 



This differs from the classical form by the addition of the exchange Pauli potential to the 
Coulomb potential, and by the modifications of the harmonic trap form in 4> c ,ext{r), Figure 
[9| If the latter are neglected it is expected that no shell structure will appear, since the 
Pauli plus Coulomb potential is still singular at short range and behaves as the classical 
mean field limit. The quantitative changes in Ciex j(r) shown in Figure [9] are modest for 
t > 2, but there is a qualitative change in shape for 2 < r* < 3 at lower temperatures. The 
corresponding density profiles from (45) are shown in Figure 12 at r s = 5. Shell formation 



in this range is clearly seen for t = 1 and 0.5. 




FIG. 12: Exchange mean field approximate density profile for r s = 5 and t = 0.5, 1, 2, 5, 10. 



IV. DISCUSSION 

There are powerful many body methods within classical equilibrium statistical mechanics 
that do not apply directly to quantum systems. Examples are molecular dynamics simula- 
tion and liquid state theory. To bridge this gap a method to define an equivalent classical 
system for the thermodynamics and structure of a given quantum system has been developed 
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PP. The objective here has been to demonstrate and test that approach with applications 
to two quite different quantum systems, the bulk uniform electron gas and charges confined 
by a harmonic trap. In the first case, a simple representation for the pair potential incor- 
porating both ideal gas exchange and RPA correlations was used within a strong coupling 
classical theory - the HNC Essential properties such as positivity are assured by the clas- 
sical formalism, in contrast to early quantum theories with mean field corrections to RPA. 
Good agreement with diffusion Monte Carlo simulation over a range of densities at t — 0. 
Results at finite t have no benchmark quantum simulation data at present for further critical 
assessment. 

The good agreement at t — is more than might be expected from the RPA input 
for the classical pair potential. In contrast, the PDW model also has similar agreement 
but is parameterized by exchange correlation simulation data at t = so agreement is less 
surprising. One possible explanation for the results here is the preservation of the exact 
perfect screening sum rule. This assures that the effective potential in the approximation 
used here has the exact Coulomb tail for large r, 

r e ( r .,^ ^ coth 2 ( ^ /2) r, (46) 

where T = (3q 2 /r is the classical Coulomb coupling constant. In particular F e (r s ,t = 0) = 

1 /2 

2T/f3hu p oc r s . It is instructive to construct a model similar to that of PDW but which 
incorporates this asymptotic property 

0(r*)^(0(r*)) (o) + A D (t,r s ,r*). (47) 

where A^> (t,r a ,x) has the Deutsch potential form 



A D («,r.,,) = ^-.*P<-3£Oj, A* =(^1 Wn («) 



Now, imposing the limit (46) gives 



Pc 



(49) 



P fihup coth (Phw p /2) ' 
Interestingly, this also determines the classical temperature. The pair correlation function 



at t — for r s = 1, 5, 10 is shown in Figure 13 



The second application here was to charges confined by a harmonic trap. The agreement 
of the model with simulation is quite good for r* > 1 is quite good, suggesting that the 
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FIG. 13: Radial distribution function g(r*) for r s = 1, 5 and 10 at t = 0. 

perfect screening sum rule captures the relevant quantum effects at strong coupling. This is 
an interesting test system for strong correlations since classically this is reflected in the for- 
mation of radial shell structure. Here, only the mean field theory (no Coulomb correlations) 
was considered as a means to explore the possibility of purely quantum mechanisims for shell 
structure. It was found that diffraction effects, which regularize the Coulomb potential at 
r = mimic classical Coulomb correlations and does lead to shell stucture. In addition, in 
the absence of both diffraction effects and Coulomb correlations, the changes in the effective 
confining potential due to exchange effects also can lead to shell structure. Elsewhere, a 
detailed application of this effective classical system will be described with all mechanisms 
for shell structure active. 
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Appendix A: RPA for the Uniform Electron Gas 

In this Appendix the RPA is defined and it is shown that the modified Coulomb potential 
reduces to the Kelbg potential for weak coupling and weak degeneracy. More generally, the 
exact large r dependence of this potential is evaluated. 
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The static structure factor is related to the Fourier transform of the hole function by 

S(k) = 1 + n J dre ik ' r /i(r). (Al) 

As a density fluctuation it is also related to the dynamic response function or complex 
dielectric function e(ui, k) [3]: 

S(k) = ^ / ^(l-e-^j^Ime^few), (A2) 

where V(k) = Anq 2 /k 2 is the Fourier transformed Coulomb potential. In the random phase 
approximation the dielectric function is 

e RPA (k,co) = l-V(k) X ^(u,k), (A3) 

and (k, of) is the response function for the ideal Fermi gas 

X<°> (k, U ) = lim / A " (£ ' k " kll) " U (£fcl) , ^ = — (A4) 

n v^o+ J (27r) rko + vr) + e kl - e| k _ kl | 2m 

1. Kelbg limit 

For weak coupling e _1 (A;,u;) can be expanded to quadratic order in V(k) to get 
S(k)^S {0) (k)-2-V(k) / dtuil-e-^y 1 (lm X {0) (uj,k)) (Re X (0 \uj,k)) . (A5) 



7T 

The real and imaginary parts of x^{w> k) are 



Re X (°) (k,,) = -^ 1 p I™ dxln(l + ze-* 2 ) (—- 

n\ 6 Ay 71 K J^oa \ / \V + K 



t\i oc jy Kj x 



(A6) 

Im x (0) k, w) = ^ — tt^P^t ln i < ^ ' A7 

The dimensionless variables k and v are 

k=- — i=i v = A = —— . A8 

A^ Ak \ m J 

Next consider the additional limit of weak degeneracy. This is implemented by an ex- 
pansion in z 

S {0) {k) -> l + 0{z), n\ 3 ^ (2s + l)z, (A9) 
Re X (0) (k, u) A + k) - ^(i/ - «)) (A10) 
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Im X (0) (k, u) -> pj^- [e^ u+K)2 - e"^ 2 ) . (All) 

with 

g(y) = -V dxe~ x2 = -2e~ u2 dxe* 2 . (A12) 

v 71 " J-oo y — x Jo 



The RPA structure factor becomes 

1 /*°° 

5 (fc) -> 1 + pV(k) = / dve~ {v - K)2 tgtv + «) - 0(1/ - «)) 

= 1 + ^(A0-^=</(V2k) (A13) 
The modified Coulomb potential at weak coupling and weak degeneracy becomes 
A ( r ) = - [ -*, e -** 



n 



-)■ - 



(2tt) 3 V^ RPA (^) ^ (0) (^) 



n./ (2tt) c 

pV K (r) (A14) 



This is the Kelbg potential of eq. (19). 



2. Large r limit 

The large r behavior of A (r) is governed by the small k behavior of S RPA (k) 

S R PA{k) ^ C oth(^). (A15) 

2mw p 2 

This is the exact perfect screening behavior [12] which is preserved by the RPA. Since S^°'(0) 
is finite at finite t and vanishes as k for t = 0, 

A(r) = — [ * e-*' / 1 ' 



r? 



(2tt) 3 V^ rpa (A;) 5(0) (ifc) 



rw P /" rfk ik . r 47r 



27™ncoth(^)7 (2tt) 3 fc 2 

r e (t, r s )r*- x , r e (t, r s ) = — \ - V, (A16) 



[3hujp coth( 



where T = (3q 2 /r is the classical Coulomb coupling constant. This is the result of eq. (13). 



[1] J. Dufty and S. Dutta, Classical Representation of a Quantum System at Equilibrium: Theory, 
previous paper in this journal. 

22 



[2] D. Kremp, M. Schlanges, W. Kraeft, Quantum Statistics of Nonideal Plasmas, (Springer- 
Verlag, Berlin, 2005). 

[3] G. Giuliani and G. Vignale, Quantum Theory of the Electron Liquid, (Cambridge U. Press, 
Cambridge, 2005) 

[4] J. Wrighton, J. W. Dufty, H. Kahlert, and M. Bonitz, Phys. Rev. E 80, 066405 (2009). 

[5] C. Jones and M. Murillo, High Energy Density Physics 3, 379 (2007); F. Graziani et al, High 

Energy Density Physics 8,105 (2012). 
[6] M. W. C. Dharma-wardana and F. Perrot, Phys. Rev. Lett. 84, 959 (2000); F. Perrot and M. 

W. C. Dharma-wardana, Phys. Rev. B 62, 16536 (2000). 
[7] M. W. C. Dharma-wardana, Int. J. Quant. Chem. 112, 53 (2012). 

[8] J-P Hansen and I. MacDonald, Theory of Simple Liquids, (Academic Press, San Diego, CA, 
1990). 

[9] K. Singwi, A. Sjolander, M. Tosi, and R. Land, Phys. Rev. Bl, 1044 (1970). 
[10] S. Tanaka and S. Ichimaru, J. Phys. Soc. Japan 55, 2278 (1986). 
[11] G. Ortiz and P. Ballone, Phys. Rev. B 50, 1391 (1994). 

[12] D. Pines and Ph. Nozieres, The Theory of Quantum Liquids, (Benjamin, NY, 1966); D. Bry- 

dges and Ph. Martin, J. Stat. Phys. 96, 1163 (1999). 
[13] C. Deutsch, Phys. Lett. A 60, 317 (1977); H. Minoo, M. Gombert, and C. Deutsch, Phys. 

Rev. A 23, 924 (1981). 

[14] G. Kelbg, Ann. Phys. 12, 219 (1963); A. Filinov, V. Golubnychiy, M. Bonitz, W. Ebeling, 

and J. Dufty, Phys. Rev. E 70, 046411 (2004). 
[15] J. Kimball, Phys. Rev. A 7, 1648 (1973); A. Rajagopal, J. Kimball, and M. Banerjee, Phys. 

Rev. B 18, 2339 (1978). 
[16] M. Baus and J.P. Hansen, Phys. Rep. 59 (1980) 1; M. Baus, Physica 79A, 377 (1975). 
[17] K-C Ng, J. Chem. Phys. 61, 2680 (1974). 

[18] F. Perrot and M. W. C. Dharma-wardana, Phys. Rev. A 30, 2169 (1984). 
[19] J. Wrighton, H. Kahlert, T. Ott, P. Ludwig, H. Thomsen, J. Dufty, and M. Bonitz, Contrib. 
Plasma Phys. 52, 45 (2012); arXiv: 1110.2465. 



23 



